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ANNOTATION 


The obiective of this study was to determine the influence of the 
channel qeometry on pressure losses and heat transfer m non-circular 
chtnnJls with hydraulically smooth walls. As a basic assumption for 
SS dJsIrlption ol “Sis influence, new integral geometrical criteria, 
selected according to experimental experience, were in ro * 

?heSe new geometrical criteria, a large set of 

pressure losses and heat transfer in circular and annular n 

with longitudinal fins was evaluated. (This research work hadbeen ? 
performed in the Institute of Nuclear Research during t y 
through 1970.) 

.„*S 

tal results of similar studies done for all kinds of cross sections. 
Detailed analysis and elaboration of experimental material men-' 

non-isothermality influence m noncircular channels. 
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1 . FOREWORD 


with nSnci?c?Lr e ??ofs n se?rinm fl? S »l, and h ? at ? ransfer “ channels 
this day oniy on e^nsl'on 

eral X ?2Ss e Sn a th° n h Compatatpon methods worked out over the past sev- 
” o St Zt r aS1 * ° f rudimenta l experimental findings a?e still 
a^Sfr- P ? by far to P rovid e with a known decree of accuracy an 

3KS to ocou? U !n ry °5 h What 1 P f? SSU P losses and fieat ?ranSfer be 

nprimLfni . cur in a channel that had not been subjected to prior ex 

US investigation . The reason for this is Constituted by the 

Sons U o fi h di »a s h®c d in these methods as basS assLp- 

SSSenSil'streSS 9 ^ 3 ^ 0 ? ° f velocl ty profiles, distribution of local 

their ICnCrIliza?ionS;C haV ? a " arrowly limited range of validity and 

nell'cS^^SSy-SiiSSSS resul.tsf ferent ° £ dhaa - 

!pec?f ic n c ;annei S gf o^ZrxteM*" 

lationslup s which describe the hydrodynamic and thermokmetic state 
m a given geometry still form an isolated point m the conStantlv el 
panding experimental materiel. Attempts to discover their mutua^rel” 
ationship have been unsuccessful. Whit we are af£Ir is nSt a compH^ 

ploliTtl 1 ilii & ii Plana ^ XOn, hUt 3 speciflc empirical method for ap- 
sible to unf?o g ? experimental materiel which would make it pos- 
e to unite all extant materiel and discover cnterial 

° £ “° re ge ? eral valldlty - in addition to Sr grCa? p?ac”CK ?mpoC P 

amSSni If °scIS«ed a and e Ci«° nShlPS '• ln diver 9 enca to the enormous 
aama n ; ot scattered and diffuse experimental materiel — would at the 

ed?«£i?S represent a basis of much better quality for studies orient- 

ed toward devising a theory of turbulent flow m channels with random 

cross sections and hydraulically smooth walls. 

However, it ought to be pointed out that the by now relatively ex 

the S rIIu?S e S mental materiel has one basic shortcoming. Comp^ison oi 
the results of measurements made by various authors on channels having 
an rdantical geometry shows in many cases that those eva^ 9 

tnr-ll en ? S that ar ? far from negligible and, at times, basic Sfflr- 

llltiin Ihllh^fS® “ e C ? n “ 3e a channel with a smooth annular cross 
section, which is the simplest geometry after a tube. So far almost 

every author who investigated the state m annular rings specifies a 

feJs S from that^f^t^ coefficient of friction losses that dif- 

s from that of other authors. A comparison of the results arrived 

K ar ^ S r!Y th °f S 13 shown ln fi g* 4 which was take-over fro! a 
tudy by Ort [4] and supplemented by results arrived at by Rybacek r 3 1 
and by the Institute of Nuclear Research. (The refatJve coe?fJc?eni 3 
of friction losses m the diagram is related to the coefficient of 
fri f t;L S? losses m a circular channel with the same Reynold's cnteri 

nSll T a S a S i ta f t i° n 1S eVSn worse ln thermokinetics of Annular chan- 
, * A detailed comparison of the results arrived at by studies of 

The t di : f? nSfer ln , thl3 JyP e of channel was compiled by Quirrenbach [191 

tens of f pe?cen?. individual Judies amoun? up to several 
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j There Can be a number of causes for such differences, and an ob- 

i! nornS^b?r nt ^ f ; act ° rs that affected individual results 
ditinnQ P anr> mo n, W J tbout a detailed knowledge of the experimental con- 
ditions and methodology used in processing of experimental data, etc. 

, it is specifically in regard to these points that most publica- 
tions provide insufficiently detailed data. puonca 

As the first attempt at generalization of experimental data can 

lnbrodu ^ tion of hydraulic diameter as the characteristic 
f f channel together with the assumption that m defining 

channel ° f % ' crit erial relations for a circular 

S? ly 1 ? bo noncirc ular channels. However, it is clear by 
. . a ^ this method has a very limited area of application and that 

-s indiscriminate use could lead to utterly wrong conclusions. At 
the same time it should be emphasized that the area of applicability 
of d h has not been clearly defined. 

The presented study represents another such attempt at general- 
ization of experimental findings. It is based on a wealth of materi- 
el obtained in the Institute of Nuclear Research during investigation 
of pressure losses and heat transfer in channels with a smooth an- 
nular cross section, annular cross section having a core with longi- 
tudinal fins and a tube with internal longitudinal fins. Approxi- 
mately sixty channels with strongly varying geometry [1], [2] have 
been investigated in the course of these experiments. Detailed know- 
edge of the conditions under which the experimental findings were ob- 
tained, together with the amount of such findings, facilitated a de- 
tailed analysis of the results and assessment of individual factors 
that could have influenced those results. In this manner were laid 
the basic prerequisites for such an experiment. 

The result of the study was discovery of integral geometrical 
characteristics which unambiguously define channels with a noncircular 
cross section with hydraulically smooth walls from the viewpoint of 
pressure losses and heat transfer. By means of these geometrical char- 
acteristics it became possible to process the experimental materiel in- 
to universally valid criterial relations for the coefficient of fnc- 
tion losses and Nusselt's criterion# 

The term "universally valid" used for describing these relation- 
ships is to be interpreted as not confining their validity only to a 
certain geometry, as comparison of these relations with published ex- 
perimental values has shown their suitability for use with channels of 
the most varied geometries that find application m technical practice. 

The author considers it to be a delightful obligation and an ex- 
pression of thanks for valuable cooperation to emphasize the signifi- 
cant contribution made by his collaborators to the presented study. 

Engr Jaroslav Hejna closed collaborated on this entire study and 
contributed a number of valuable suggestions and critical comments to 
its final shape. 
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2. LIST OF SYMBOLS 


X 

longitudinal coordinate 

[ m ] 

ll 

height of fin 

[ m] 

i> 

width of fin; rod spacxng 

[ m ] 

n 

number od fins; exponent 

[1] 

d 

diameter of circular channel, rod diameter 

[m] 

V 

root diameter of finned tube 

[m] 

<M 

i—i 

inner, outer annular diameter 

[m] 

l 

length of model; length of measured section 

[m] 

Ij 

length of heated section 

[in] 

l H 

length of hydrodynamic 

[m] 

l t 

length of thermal 

[m] 

S 

wetted circuit 

[mj 

F 

cross-sectional area of flow 

[m2] 


fin circumference 

[m] 

d k 

hydraulic diameter 

[mj 

? 

specific weight of gas 

[ kg/m “ 3 ] 

X 

thermal conductivity of gas 

I W/m“ Vdeg - 1 ] 

Xm 

thermal conductivity of model tube material 

[W/m~Vdeg -1 ] 

]X V 

dynamic, kinematic viscosity 

[N/s/m“ 2 ; m 2 /s 

C P 

specific heat at constant pressure 

[ J/kg -1 /deg _1 ] 

P 

static pressure 

[ N/m“ 2 ] 

t 

local temperature; temp, at root circle 

[ deg] 

0 

median gas temperature 

[deg] 

* 

thermal factor 

[1] 

u 

velocity 

[m/s -1 ] 

u + 

dimensionless velocity 

[1] 

y + 

dimensionless distance from wall 

[1] 

U 

friction velocity 

[ m/s -1 ] 

G 

flow-through weight 

[ kg/s" 1 ] 

M 

overall thermal output 

[w] 


thermal output per unit of length 

[ W/m -1 ] 

a 

coefficient of heat transfer 

[ W/m -2 /deg -1 ] 

X 

tangential stress on wall 

I N/m" 2 ] 

n? 

fin efficiency 

[1] 

L 

see relation in definition (29) 

[mj 
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L 


see definition m equation (29) 

[m] 

S * I 


see definition in equation (31) 

[m] 

1 = r 


integral geometrical criterion 

tl] 

c r° 

^A> * H> 


see para 6.2 


H 


geometrical factor 

[1] 

e 


coefficient of friction losses 

(1] 

Mu, Re, 

Pa. 

criteria of Nusselt, Reynolds and Prandtl 

[1] 

M 


median value of factor N along 



cross-sectional area of flow 
or along 

W mux maximum value of factor N in channel 
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IN NONCXHCuLTS^f wM^ l “LoOTH WALLS 

J* Malak 

Instituted Nuclear Rfsearch 063 ! 5 rague ' Czechoslovakia, 

ear Kesearch - Information Center in Rek 1 

=ie« r RlIea^r&1 s ^ P L ri ^ Institute of ■_ j 

losses and^eat transfe^in oh exp f rimen tal research into pressure 

having longitudinal filsdhichd? Cir ^ ular and annular cross 

lndm ea d eS6arch in the years'l967l70 L IS at the In stitute 

ciiicl [ 2 ]« These reports ai en t • . Provided in reoo:rt «5 rii 

form of relationships 1 1St the resuits of measurements in thl 


5 = i (Re) and Wu = ^ 


(Re) 


in £sris$£; ti £ ssinSM Si s u st ° u f d - y . 

a rectangular cross^sectior^with l JttTn7 Vo The flns ^ere of 

was a total of 27 variants which leJSfsentlT V OTT1 * Investigated 
values of three basic geometrical la SIlISl combinations of three 

lected: channel width (or cCCCI sill? ? S r Which they were se- 

fms and height of fins. The length of^la ° U ? er tube D 2 ), number of 
cases and was equal to 1,200 mm Inllt fntTtt W * S ldent ^al in all 

in 1 fig?°?? d affect the results of measurements by^l^hlpe^sdSSI 1 

inner logitudinal fins^Vrhis oeomet™^^!? 11111 Clrcular channels with 
B.) The fins were eSintiSu? ZhTstnl referred to as geodry 

larly to geometry A, it was nlann^a ? mS 3S those ln geometry A. Simi- 
mg to combinations of three values n! f xamin e 27 variants correspond- 
ed height of fins. HoleleVliree vaATts 2t“ te f ° f tube ' »p5c±Sg 
all, as finned tubes for those varianfc 1 ^^ 616 not investigated at 

P | eS | ribed tolerances and with acceDtablldd^ be P roduce d with the 
of channels in all variants was I surface quality. The length 

inlet is shown xn fig? 1 ?? S WaS 2 '° 50 - nUn * c °nfiguration of the chllnel 


li ra t^ O ^ aaSe ^ ea ^ ea ^uIsr e channels P ?i e 200 te<3 1^ 1 r GStlgati ° n ° f p — 

limiting cases of geometry a at h-*0. ' ° i° n 5) which represented 
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Detailed data regarding all of the investigated channel models 
are listed in appended tables. 

All experiments were conducted in an open air line with a venti- 
lator located behind the outlet from the experimental section. The 
models were heated by resistance heating at a constant thermal output 
per unit of length of the model . 

The measurements included m all cases the flow-through weight of 
air, thermal output, distribution of static pressure along the channel, 
distribution of temperatures of the heat-exchange surface on the root 
circle along the length of the heated model and air temperature at the 
input and output of the experimental section. 

Investigation of pressure losses was done m all cases at isother- 
mal flow. The friction loss coefficient was computed with respecting 
the difference of the pulse at input and output of the measured sector 
according to the equation 

5 = 2 d h !f - i] (1) 

wherein _Aj3 = pressure loss in channel sector of length 1 

p,£ = mjsdian^values of pressure and specific weight of air m 

G = flow-through weight 
F = cross-sectional area of flow 


Reynold's criterion is given by the relation 

Re = f - ‘ - - d K (2) 

For computation of the median value of the coefficient of heat 
transfer a (centered along the heated circumference) was used equation 

qi = a(ir.t> - n.4 + - 9) (3) 

Derivation and analysis of the validity conditions for the equation/9 » 
is contained in [1]. Of importance to us is the fact that in our exper- 
imental conditions this equation is valid irregardless of the actual 
distribution a along the heated surface. 

i 

The connotation of individual values in equation (3) : ^ 

n, & = number and width of fins 
V = root diameter of fins 
S* = fm circumference 
Z = wall temperature at root circle 

0 = median temperature of gas 1 


Efficiency of fins with a rectangular cross section was taken ac- 
cording to Gardner [16]: __ 

tgh w , , ^ 2a 

where X m is thermal conductivity of the fins' material. 


(4) 


In determination of the delivered thermal output from unit length 
of the model q ,, m view of the structure of the used models (thick- 
walled tube with fins) consideration had to be given to axial heat son- 
ductance. In view of the high thermal conductivity of the material of 
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the finned tubes and the structure of the heating system we can as- 
sume an automodel distribution of temperatures m radial cross cuts 
of the model. This can be put m the form of 

_ N . e . 2 z t , c . 

*L + Fm * X ' JI 3P r (5) 

The median temperature of gas at distance from input into the 
heated model at random distribution is determined by the relation 

0 = 0o + / A- d x (6) 

yj • c-p 

where W = overall thermal output ,. 0 

/y = length of the heated model ' 

F m = cross section of finned tube 

0 O = gas temperature at input into experimental section. 

A program was compiled for analysis of experimental values on 
the computer GIER. To facilitate solution of the system of equations 
(3) t (6) , the experimentally determined distribution of temperatures 
along the length of the model was substituted for each regime by the 
method of least squares with a polynomial of the fifth order. At the 
same time, for the sake of comparison, was also determined the median 
coefficient of heat transfer without taking into consideration the ef- 
fects of axial heat conductance. Results of thermokinetic research 
were processed in the form of the relation 

Hu • i (Re) 

where values of Reynold's and Nusset's criteria are counted with values 
of thermophysical properties of the coolant in relation to median tem- 
perature and gas pressure in the given radial cross cut of the experi- 
mental channel. 


In addition to 3T, Nu and Re, in all cross cuts with measured sur- 
face temperature was also determined the thermal factor ip 

* * § < 7 ) 

4. Effects Distorting Experimental Results /I 

The basic prerequisite for succeeding m out attempt to unify and 
generalize the experimental material was to base this study only on 
cases in which measurements were not subjected to distorting effects. 
Talking of factors that can distort experimental results, our natural 
premise was that the objective of our experiments was to measure pres- 
sure losses or heat transfer in the region with developed hydrodynamic 
conditions (distribution of velocity, distribution of tangential stress 
on channel walls) or, eventually, with developed thermal conditions (au- 
tomodel thermal profile (£-0) ) in channels with hydraulically smooth 
walls . 


Distorting factors can essentially be divided into two categories; 
The first category includes factors which m their essence are bound 
with the examined physical process and their distorting effect must be 
viewed as the consequence of our insufficient knowledge of the examined 
process. They are: 

1. transitory region of flow 4. non-isothermity of flow 

2. hydrodynamic entry-end sector 5. roughness of walls. 

3. thermal entry-end sector n 



In addition to these factors , in individual cases come more-or- 
less into play factors of the second category which are a consequence 
of the selected experimental concept, failure to maintain the speci- 
fied boundary conditions in its implementation, lack of suitability 
or inaccuracy of the used measuring methods, etc. Assessment of such 
distorting effects (which should rather be called experimental errors) 
and an eventual making of suitable corrections m processing of exper- 
imental data is, however, within the realm of possibility of only and 
solely of the author of any given experiment. 

Let us have a closer look at the factors in the first category, 
which in our opinion are the key reason for a number of contradictory 
data found in experimental materiel. We attempted to cull from all 
accesible references the current state of knowledge regarding these 
problems and supplement them by our findings as well as by the re- 
sults of an analysis of our own experimental materiel and that of 
others. Our objective was to redefine more accurately individual 
terms as well as to qualitatively assess the effects of individual /12 
factors in the case of noncircular channels. 


4.1 Transitory Region of Plow 

The term "transitory region of flow" is used m hydrodynamics in 
two different connotations. In most cases we associate this term 
with that part of dependence of the coefficient of friction losses 
in a circular channel on the Reynold's criterion, m which £ m tran- 
sition from dependence on laminar flow to dependence on turbulent 
flow increases \*ith increasing Re. For the sake of differentiation 
we shall refer to this region and the transitory region 


We encounter this term m its second connotation in describing 
the velocity distribution m turbulent flow. This distribution is 
usually expressed m the form of the relation 


u 


+ _ 


- i iy + ) 


u 


+ _ 


a 

77 *- 


.* - 


/t 


y 


y.v' 

V 


jr = friction velocity 


( 8 ) 


Here u is local velocity ay distance y from channel wall and x is 
tangential stress on the wall. 


In the area of the laminar substratum applies that 

u + = y + (9) 

In turbulent core of flow in a circular channel — and, accord- 
ing to a number of experimental results also in noncircular channels 
— the distribution of velocity occurs m accordance with the relation 

u + = a. tog y + + b (10) 

where constants a, b according to some authors have universal validity 
while, according to others, they depend of the geometry of a channel. 
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Between these two regions lies the region of transient flow — /13 

which we shall refer to as transient region a + — m which occurs a 
laminar as well as a turbulent transfer ofthe momentum of motion. 
However, according to the latest findings this process cannot be in- 
terpreted as simultaneous effects of these two types of transfer, but 
as an irregu f ar alternation of laminar and turbulent flow. A factor 
called coefficient of intermittence (discontinuity) y — which is de- 
fined as the time component during which turbulent flow prevails at a 
monitored point of the examined space — is used due to its vividness 
for description of transient flow. Thus, if applies that y = 0, the 

flow is purely laminar, while at y = 1 the flow is completely turbu- 
lent. 


In a circular channel the extent of the transitory region £ (ex- 
pressed by the extent of Re) as well as its lower boundary depends on 
a number of factors, the most important among them being geometrical 
conditions and the degree of turbulance at entry into the channel. 

In routine practice it is assumed that the transitory region lies at 
Re = (2,3 + 3, 5). 10 . Nevertheless, during special conditions at en- 
try, it became possible to maintain laminary flow into substantially 
higher values of Re. 

The best insight into the processes that occur m the transitory 
region B, in a circular channel from the viewpoint of velocity condi- 
tions is provided by the results obtained by J. Rotta [15]. From de- 
measurements of velocity states and their chronological changes 
he found out that in the transitory region B, throughout the flowthrough 
cross section the velocity distribution irregularly transfers from a 
developed laminar profile into a developed turbulent profile and vice 
Hereby the coefficient of intermittence is not dependent only 
on Re, but simultaneously undergoes strong changes with distance from 
entry into the channel (fig. 5). Thus, in a cicular channel the tran- 
sitory region £ corresponds to velocity regimes during which the tran- 
sient region u+ affects practically the entire flow-through cross sec- 
tion of the channel. 


Velocity conditions in a noncircular channel will differ from 
those in a circular channel in direct relation to the complexity of 
channel configuration. As a measure of complexity of configuration we 
can use the ratio of the maximum and median width of channel in the 
direction of normals to the surface, or the ratio of the maximum and 
median curvature of the wetted circuit. The most suitable and most 
objective measure of the complex configuration of a channel is distri- 
bution of tangential stress on the wall, or (T max /t) . 

In channels with complex configuration we do not deal any more with /I 
a cross section represented by a single simple region as is the case in 
a tube. In extreme cases we could rather talk of several regions that 
strongly differ from the hydrodynamic viewpoint and that retain a cer- 
tain measure of "autonomy". For better visualization m limiting cases 
we can think of individual regions as independent channels connected in 

The hydrodynamic conditions in each of them would then be 
defined by its own Re. 
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In a tube during developed turbulent regimes the laminar sub- 
stratum region as well as the transient region u + affects only a 
negligible part of the flow-through cross section. However, m chan- 
nels with a highly complex configuration, laminar flow can be retain— 

111 some parts of the flow-through cross section up to relatively 
high values of Re, irregardless of the fact that in the remaining 
parts of the cross section prevails parallel turbulent flow, whereby 
the share of regions with laminar flow can be relatively high in re- 
gards to the overall flow-through cross section. 

These regimes are vividly reflected in the results of measure- 
ments done by Eckert and Irvine [30] who investigated conditions in 
a channel with a cross section in the shape of an isosceles triangle 
with an apex angle of 11.5°. By means of visualization of flow and 
measuring of velocity profiles they delineated m the cross section 
the area of laminar flow and the dependence of its extent on Re (see 
fig. 51). At a still relatively high value of Re = 10 4 , at the apex 
of the triangle, up to yg- of its height prevails purely laminar flow. 

The compounded structure of flow in a single channel, or simul- 
taneous effects of laminar and turbulent transfer of motion necessar- 
ily must become reflected also in a loss of pressure. Under these 
conditions the dependence of the coefficient of friction losses on Re 
will be analogous in character to that of £ in the transition area m 
a tube. Thus, from the viewpoint of local velocity conditions, the 
term "transitory region £" assumes a broader connotation in the case 
of channels with complex configuration than was the case with a tube. 
However, from the viewpoint of physical substance the connotation of 
this term is the same for all channel geometries. Thus, under the 
term transitory region £ we shall subsequently understand that part 
of Re in which m a given channel the friction pressure loss is the 
result of the effects of laminar and turbulent transfer of motion. 

On the basis of what was said about velocity conditions m non— /15 
circular channels it could be expected that the transitory region in 
noncircular channels will be wider than m a circular channel and that 
the lower boundary of the region of purely turbulent flow will be 
shifted toward higher values of Re. This assumption is fully confirm- 
ed by experimental results . 

From the hydrod-namic viewpoint, the noncircular channel closest 
to a circular channel is an annular channel. The results of our mea- 
surements of pressure losses m six annular channels with a curvature 
of Dr/D 2 = 0.65 to 0.95 show that great differences exist already be- 
tween circular and annular channels. Our results are plotted in fig. 

3 in the form 1/I£ = & (Re^) supplemented by theoretical dependence of 
laminar flow in annular channels with V\ /V 2 = 1. For the sake of com- 
parison the diagram also shows the corresponding dependences for a 
smooth tube. It is obvious that the transitory region in the case of 
an annular channel is substantially wider than that of a tube. When 
we attempt to estimate from the progression of the dependence//^ = 

/ [Rt Cl, interspaced by exp. points, the limits of the transitory 
region, we obtain the values 

Re s 2 . 10 3 
Re"> 2.4 . 10 4 
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lv 6 ^ 111 ^ 85 ° and H denote the lower and upper boundary, respective- 
vY' . In ^termination of the upper boundary it is questionable whether 
£« a nS !- nt ® rs P e f sed Wlth experimental points does really correspond 

pure y urbulent regime, or whether at higher values of Re the 

?h^ a o?f° g J eS ^ mi ?? t not resemble more closely, e.g., the relation- 
snip uij which will be discussed m detail sub para 6 and 7. The 

± a ct that in annular channels the boundary between the transient and 
turbulent region is less conspicuously defined than in the case of a 
tube, is the logical consequence of the already mentioned difference 
velocity conditions. Nevertheless, it is clear that the region of 
purely turbulent flow m annular channels lies roughly by 1 order of 
Re higher than in the case of a circular channel. 

Results analogous to those for annular channels were also obtain- 
ed in our measurements on channels of geometry A and B which represent 
channels with one of the most complex configurations that have been 
investigated so far. Among our results there are a number of cases 
with a significant share of the transitory region, or even completely 
dominated by the transitory region. Figures 11 and 12 show several /16 
such typical cases for geometry A and B. 

The effects of geometrical conditions at the entry into the chan- 
nel on the extent and boundaries of the transitory region £ can come 
to the fore very strongly even m channels with a complex confiqura- 
tion. Figures 7 and 8 show the plotted results of our measurements - 
of four variants of geometry A at varying entry conditions. Each of 
the measurements was done using two models, in spite of the fact 
that these models were to be absolutely identical from the hydrody- 
namical viewpoint, measurements of pressure losses differed. Regret- 
ably, when we found out about these differences, the models II had 
already been destroyed and we could not undertake a detailed compar- 
ison of both models. However, supplementary verification of the pro- 
duction process turned up the fact that in some cases the edges of 
the model at its entry were rounded while m other cases their sharp 
edges were retained. Experiments with models I were therefore repeat- 
ed under modified conditions at their entry, which we achieved by 
placing sheet-metal screens m front of the entry as is shown m fig. 

1. In this manner we obtained a third group of results which, in some 
cases approximate more the original results obtained with model I, m 
other cases the results obtained with model II . For enhanced orien- 
tation and easier comparison the results have again been plotted m 
fig. 9 and 10 in the form 7/I£ = i (Rel£). From the progress of indi- 
vidual measurements it can be deduced that with increasing value of 
Re these dependences transform into a single dependence. 

Analysis of the results obtained by other authors in measurement 
of pressure losses in channels of the most varied geometries led us 
to the conclusion that a number of those results were obtained in the 
transitory region ? without the authors of those studies being aware 
° f lu 1S f consequence of incorrect application of findings re- 

garding the extent and boundaries of the transitory reqion l in cir- 
cular channels to noncircular channels. That, as a result, leads to 
a number of incorrect conclusions deduced from experimental results 
affected by the transitory region. 
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•f n conc ^- us;LOn we ought to mention one more factor on the effects 
of which we were withdrawing comment for the time being. The results 
obtained by J. Rotta clearly show that stabilization of flux condi- 
tions in the transitory region £ m the case of a tube calls for ab- 
normally long channels (see fig. 5 and 6) on the order of several hun- 
dred a fo. The situation will be even worse m the case of noncircular /17 
channels. As channels used in investigation of pressure losses never 
attain such lengths, the results of measurement m the transitory re- 
gion will probably be distorted even more by the effects of the entry 
sector. 


4.2 Hydrodynamic Entry Sector 

The term "hydrodynamic entry sector" denotes that part of the 
channel beyond the entry in which occurs development of the velocity 
profile from the initial profile at entry up to a fully developed 
profile corresponding to the preset hydrodynamic regime. Connected 
with the development of the velocity profile is also development of 
tangential stress or its distribution along the wetted circuit, and 
development of transfer characteristics. The length of the entry sec- 
tor depends on the hydrodynamic regime, the geometry of the channel, 
configuration of entry or, eventually, on defective flow prior to the 
experimental part. 

Let us first examine the case of a tube. For laminar flow and a 
sdged entry (i.e. , with a balanced velocity profile of input 
velocity) we theoretically derived the relation 

• Re (11) 

The constant of proportionality differs from one author to another and 
ranges between 0.03 to 0.06 (see, e.g. , [11]). 

In the case of turbulent flow, the decisive effect on the length 
of the entry sector accrues to the shape of the channel entrance. It 
is this shape that determines whether the boundary layer in the devel- 
opment of the velocity profile will be turbulent from the start, or 
whether at the beginning part of the channel it will first be laminar 
and only after a certain distance from it transforms into a turbulent 
boundary layer. As a rule of thumb, the aerodynamically smoother the 
entrance, the longer the entry sector. From experimentation we know 
that an aerodynamically smooth entrance calls for a channel length of 
least 50 d if we are to achieve a perfectly developed velocity pro- 
file. The effect of Re are for all practical reasons absent. 

As we already mentioned, the situation in the transitory region E 
is entirely different, as is also shown by the results obtained by J. 
Rotta. Figure 5 shows the dependence of the coefficient of intermit- 
tence in the transitory region on the relative distance from the en- 
trance into the channel, and fig. 6 shows change of velocity m the /18 
axis of a tube with Re for 3 distances from the entrance. Both figures 
show that stabilization of the velocity profile m the transitory re- 
gion calls for channel length on the order of hundreds of diameters 
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In case of a geometry most closely resembling a tube, i.e., a 
smooth annular channel, the situation m the turbulent region is not 
clear. On the one hand, only very few experimental data are avail- 
able and, on the other hand, these data show basic differences. In 
in-\esti gating pressure losses m annular channels m a range Vi /£> 2 = 

= 0.77 through 0.93, J. Ort [4] measured the values of the entry sec- 
tor m a range of L^/ d[ t = 7 through 21. In keeping with his results 
is also Dreissler's datum for a flat channel Ltf/d/j = 10 (quoted from 
Ort's study). 

Results analogous to Ort's were also obtained m our experiments 
with models of geometry A and B. The values of Ltf/dh ranged between 
5 through 30, whereby a majority of cases were at the lower limit of 
that range 1 . 

Extreme maximum values of L^/d^ were determined for variants in 
which the investigated range of Re was completely or partially in the 
transitory region £. 

We were assessing the length of the entry sector, similarly to 
Ort, from the progression of static pressure along the channel, or, 
to be more specific, from a comparison of the actual distribution of 
pressure with linear progression which was interspersed by experiment- 
al points in the region without apparent effects of the entry sector. 

The dijference between actual and linear progression was in most cases 
very small. In view of the aforementioned and with a view to the rel- 
atively short length of the entry sector (and the resulting small num- 
ber of experimental points in the entry region) this procedure does 
not appear to be very accurate and the quoted values of Lft/d/, must be 
regarded as only approximative information. 

However, on the other hand, Rothfus with his collaborators [41] /19 
in experimental investigation of annular channels m the range of 
V\ /Vz =0.12 through 0.35 came to the conclusion thatstabilization of 
hydrodynamic conditions calls for a channel length in excess of 200 dfo. 
They based it on direct measurement of tangential stress on the inner 
diameter and its changes with distance from the entrance into the chan- 
nel. 


From what was said about experimental techniques used in individ- 
ual cases we get an explanation of the difference between our results 
and those of Ort on the one hand, and Rothfus' results on the other. 
While Rothfus based his approach on measurement of local values, our 
and Ort's measurements provide information only about the median value 
of tangential stress on the wetted circuit. And that is a substantial 
difference, the consequences of which we will best realize by a compar- 
ison of conditions m a circular and a noncircular channel. 

*1 ■ ■ 

In keeping with these values are also the results of experiments on 

enlarged models of geometry A at Vi /V 2 = 1 that are currently conduct- 
ed at the Institute of Nuclear Research. The accuracy of all measure- 
ments is at the same time substantially higher than was the case m 
experiments described m the report [1]. 
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Let us start out with a tube. We assume the presence of balanced 
entry conditions at the entrance into the channel. Under this assump- 
tion the local tangential pressure on the wall is identical with its 
median value along the circumference of the channel. With increasing 
distance from the entrance the gradually balanced velocity profile 
transforms into a developed profile. In the region close to the wall 
the velocity decreases while, conversely, it increases in the central 
part of the channel. There is a drop in the radial gradient of veloc- 
ity at the wall and, as a result, with increasing distance from the en- 
trance there is also a monotonous decrease m tangential pressure on 
the wall. 

The situation is substantially different in the case of a noncir- 
cular channel. Let us again assume a balanced velocity profile at the 
entrance, meaning that input velocity coincides with the median veloc- 
ity in the channel. However, during development of the velocity pro- 
file the situation is not as unequivocal as was the case in a circular 
channel. In view of the configuration of the channel m some regions 
of the flow-through cross section the "median" velocity will decrease 
while the opposite will apply to other regions.. Tangential stress on 
the wall in individual regions will also change correspondingly. From 
our and Ort's experiments we can draw the conclusion that development 
of hydrodynamic conditions m channels with cross sections of a com- 
plex configuration (i.e., with x m ax/i > 7) progresses m such a manner 
so that — with the exception of a small region beyond the entrance cor- 
responding to values of Ltf/dh from our and Ort's experiments — the 
development of local velocity conditions and local valus of x has no 
effect on the median value of T and the latter, for all practical pur-/20 
poses, does not change along the channel. According to our results, 
this conclusion should essentially apply also to the transitory region 


4.3 Thermal entry sector 

The term "thermal entry sector" or, also, "sector of thermal sta- 
bilization" denotes for us that part of the heated channel in which 
occurs the development of thermal conditions from the initial balanced 
thermal profile in coolant at the entrance to the heated section up to 
the region of automodel of thermal profile (t - 0 ) . 

In implementing an experiment or the real thing there can occur 
two different cases : 

1. The beginning of the heated section is identical with the en- 
trance into the investigated channel. Development of velocity and 
thermal conditions progresses simultaneously. 

2. Ahead of the heated section is an unheated channel of ident- 
ical geometry and, as a result, at the beginning of the heated sec- 
tion the velocity profile is already fully developed in isothermal 
flow. Another relatively small change in hydrodynamic conditions oc- 
curs only in the flow ceasing to be isothermal. 
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4 z^ 0 +. AS in act ^ a ^ Practice m a majority of cases is encountered the 

' and be !f ause our experiments were also conducted under such 
tho d J^ 10 ? 8 ' We ^. Sh f 1 i aubse< ?uently deal only with this arrangement. As 
the development of hydrodynamic and thermal conditions progresses in 
the entry sector simultaneously with mutual interaction, these two 
processes cannot be separated from each other. Thus, whenever we re- 
fff t ° the thermal entry sector we will be referring to that part of 
the channel where there occurs development of flow simultaneously with 
the development of thermal conditions. 

According to experimental findings, the lengths of hydrodynamic 
thermal entry sectors m a cicular channel are for all practical 
purposes approximately the same [11]. Most authors who have experi- 
mentally investigated this problem specify values of Lr/dt* ranging 
between 40 and 60. In the case of noncircular channel!: sam! S was 
the case of local hydrodynamic conditions, substantially higher values 
can be expected to be encountered than was the case with a tube. A 
few experiments conducted in this respect fully confirm such an as- 

• _ Tbas ' e *9*» Eckert and Irvine in their already quoted study /21 
[30] specify for the investigated triangular channel L T % 100 du. The 
effects of the entry sector together with those of the transitory re- 
f lso became strongly reflected in most published studies devo- 
tad Z° beat trans fer m annular channels that used the method of con- 
stant temperature of the heat-exchange surface in determination of the 
median coefficient of heat transfer m the entire channel. Assessment 
e 5T ors and limitations of this method from the viewpoint of deter- 
mination of a in the region with developed velocity and thermal condi- 
tions was done under the assumption applicable to a tube. This would 
also account for the great differences between the results of experi- 
ments conducted by various authors (see [19]). 


We based our assessment of the effects of the thermal entry sec- 
or during our experiments with channels of geometry A and B on mea- 
surements of the progress of temperature of wall along the length of 
the model and the corresponding distribution of a. To enhance clarity 
and suppress the effects of Re the results were processed in the form 


*= * 


(x) 


(12) 


For reference value a ^ was taken the value of a in the last measured 
cross cut near the channel exit. 

In assessment of the processed results, m the case of some vari- 
ants there appeared certain errors and lack of logic m the distribu- 
t 5 ab we f e ^ aused b y ioss of heat at the end of the model. 
S thfllnlfnn ?£ r i? c J u ?i on . of the effects of axial heat conductance 
effect*? nf ? he dlstributlon of did not manage to eliminate the 

effects of heat losses m some cases. This calculation method — de- 

scribed m para 3, is based on experimentally determined distribution 
of temperature m the wall of the model, which by the method of least 

^ are f 5 eplaC K S , by a P ol y nomial bb e fifth order. This polynom- 
ial, . which described temperature distributions with a high degree of 

n0t auff;Lce ln some cases for the first and second de- 
vation of temperature at the ends of the model. In such cases the 
assessed progress of a at the end of the channel again rises. Then it 
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became necessary to carry out an analysis and comparison of progres- 
sions of IT, evaluated with respecting as well as neglecting the ef- 
fects of heat conductance, compare the progressions a/ a„ for various 
values of Re and determine to what distance from the end the effects /22 
of heat loss extend. Such an example can be seen, e.g., in fig. 14. 

? f ac ili tate comparison, fig. 13 and 15 show cases when the effect 
of terminal losses failed to materialize. 

We determined that the length of the thermal entry sector in our 
experiments ranged between L-r = [17 -SO) cf/ t according to the type of 
the examined geometrical variant. As the length of the models ranged 
between 50-400 the region with developed velocity and thermal con- 
ditions was never reached in the case of some variants. 

In assessment of the effects of the transitory region and of the 
entry sector, very effective proved comparison of the relations 

|=<(Re/?| [? ) 

It' <UI M 

Mu. = t HRe) [Mu] 

Such a comparison and placement of the relationships (?) , (a) and [Mu) 
into mutual connection provides a more comprehensive idea of the hyd- 
rodynamic and thermokinetic conditions in the channel and their change 
along the channel. Let us cite a few typical cases: 

In variant 13 the relationship (?) (fig. 11) up to a value of Re 
equal to 8.10 has the character of a transitory region. At the same 
value of Re also the relationship (a) (fig. 13) shows a strong region 
of laminar heat transfer, and the relationship [Mu) (fig. 18) , only at 
values of Re > 8. 10 3 transforms into a relation of the type 

Mu = . Re 0,8 

Analogous facts could be found m the case of variant A 15. 

In the case of variant A 25 (fig. 11, 16, 18) all three relation- 
®nips lie clearly m the transitory region and closer to the laminar 
region. (Values of Mu in fig. 18 were assessed in the one-before-last 
measured cross cut only for illustration of the conditions.) 

, Variant B 18 represents a case in which there occurred a distor- 
tion of the relationship (a) (fig. 17) as a result of a terminal loss, 
primarily m the last cross cut. For the progression of the relation /23 
Ca) at various values of Re it is readily obvious that in the region 
. < 3.10 practically the whole model lies m the entry sector. In 
view of the fact that relationship (?) (fig. 12) up to Re = 4.10 4 is 
in the transitory region, the great length of the entry sector is to 
be expected. The relationship (Mu) also behaves in keeping with 
these factors (fig. 19) . 

These several typical examples clearly show how interrelation of 
he relationships (a) , (?) and (A/a) makes it possible to identify even 

the relatively weak effects of distorting factors. 
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In conclusion, the effects of the entry sector (or effects of the 
length of a model) on the relationship Mu = (Re) can be most vividly 
demonstrated in fig. 20 which shows the results of measurements on the 
variant B7 as a typical example. The points marked by circles are lo- 
cal values of Wu at varying distances from the entrance into the chan- 
nel plotted in relation to the local value of Re. The almost vertical 
intermittent lines represent the progression of Mu aLong the length of 
the model at constant flow of coolant. Values of Mu decrease with in- 
creasing distance from the entry and at balanced conditions gradually 
attain a value corresponding to the sought after relation Mu = (Re) . 
The second system of intermittent lines shows how would progress if 
we used models of varying length m our experiments. The resultant 
relation Mu = j((Re) is then distorted by the effects of the entry sec- 
tor not only in regards to absolute values, but also form. The result- 
ant form of the relation Mu = $ (Re) for the variant B7, shown by the 
thick solid line throughout the investigated extent of Re, corresponds 
to the typical form of this relationship for a circular channel with 
inclusion of the transitory region of flow (see, e.g., [3]). 


4.4 Effects of Nonisothermal Flow 

In assessment of thermokmetic experiments we relate all thermo- 
physical properties of the coolant to its median temperature, i.e., 
temperature computed from thermal balance. It appears, however, that 
with increasing difference between the temperature of the wall and the 
median temperature of the coolant (i.e., with increasing thermal flow) 
at otherwise identical conditions of the experiment, the value of the 
thermal coefficient decreases. This phenomenon is caused by a change 
in the thermophysical properties of the coolant m the region at the /24 
heated wall as a result of high temperature gradients. This means 
that median temperature of coolant is not suitable as the defining 
temperature for thermophysical properties and that the correct value 
is to be found somewhere in the region between t and 0. This fact is 
respected in the critenal equation for Nusselt's criterion by a suit- 
able function of the so-called temperature factor ip which is defined 
as the ratio of wall temperature and the medium temperature of coolant 
(see equation (7)). 


Kutateladze 12 derived on the basis of an analogy between trans- 
fst of motion and heat for gas— cooled circular channels the relation 
between the coefficients of heat transfer at isothermal and noniso- 
thermal flow in the form 


a , 2 f 

a 0 'Vty + r 

wherein a 0 is the coefficient at =1. 

However, most empirical findings are expressed in the form 


(13) 

(14) 


Relation (13) can then be replaced by relation (14) with n = 0.55. The 
values of n from the experiments range between 0.3 to 0.55. However, 
we must realize that conducting an experiment that would facilitate di- 
rect assessment of the effects of the thermal factor is for all practi- 
cal purposes impossible. In actual implementation of an experiment of 
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this kind a change m is immediately accompanied by a change 

and Re. As a result, assessment of the effects of is always m 1 - 

rect and burdened with large errors. 

Since we are unable to describe the effects of n ° nia °“ e Jf“r“?d 
with any degree of accuracy even in a circular channel, then in regard 
to channels with a complex configuration, where temperature of t 
heat-exchange surface varies along the circumference, we know P r a°tl 
cally nothing about these effects. At the same time, the error. cause 
by the thermal factor can be considerably large. In the cas 
tube? e.g., change in * from a value of 1 to 1.3 corresponds to a 
change in Nu by approximately 13%. 

For that reason, operational conditions investigation of heat 
transfer in Geometry A and B were selected so that the thermal racror 
in the region of balanced thermal conditions be as identical as Po " / 
for all reaimes of all tested geometrical variants. The thermal 
?a“or in cur experiments is related to the wall temperature on the 

root circle. 

We also attempted at the same time (as we shall mention later) to 
experimentall^determine this effect for at least several variants. 


4.5 Effects of Wall Roughness 

ano-fher factor that can strongly affect measurement of pressure 
«» .-.c The effects of roughness are relative- 

meter k/P, (fe = height of unevennesses) . We know substantial y 
in regards to cases when tte density of unevenesses is low and when 
the effects are produced by their geometrical con igura ion. 

in variegated channels the same absolute roughness will manifest xt- 
selFJn Sd?vidua? parts of the channel m a different manner. 

We only know that, same as in case of a tube, manifestation^ 
the effects of wall roughness increases with increasing Reynola - 
criterion? As soon ns wall roughness starts to affect the flow - 
ditions the gradient of the relation Z - 5 (Re) starts to na.oi a 
in the end, this relation will transfer into the automodel region. 

Rendering an objective decision whether the effects of roughness 
became manf felted in a specific case - with the exception of the 
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au ^' orno< ^ e ^- region is impossible as long as we do not posses a de- 
tailed knowledge of circulation in noncircular channels. At the pre- 
sent time such assessment depends on the subjective explanation of 
the experiment 1 s author or, more precisely, on the criteria selected 
by him for the assessment. 

According to the criterion selected by us, no effects of rough- 
ness were encountered in our experiments. 


5. Basic Findings and Assumptions 


/26 


5.1 Pressure Losses 

From a detailed processing and comparison of experimental data 
from our own as well as others' experiments, accompanied by a detail- 
ed analysis of experimental conditions, we reached the conclusion 

in all cases, in which there is no manifestation of the formerly 
mentioned distorting effects, the coefficient of friction losses in 
the region of developed turbulent flow corresponds to the general e- 
quation 

71? = Ci .log Re/? + C 2 (15) 

wherein constants C\ and C 2 are the function of channel geometry. 

The empirically proven general validity of equation (15) thus 
narrows down the problem of finding a criterial relationship for the 
coefficient of friction losses to determination of the dependence of 
constants Cj and C 2 on channel geometry. However, from this view- 
port' processing of experimental data in the form of the relation 
'//? = 6 (Re/5) proved unsuitable. In this system the effects of ex- 
perimental errors is reflected much stronger than in the system ( 5 ; 

Re ) , resulting in a larger scattering of experimental points. In 
view of the fact that our experiments in individual cases did not cov- 
er even one whole order of Re, we found it very difficult and problem- 
atic to accurately determine constants C a and C 2 , which are very sen- 
sitive to interpolation of the relation (15) . This applies primarily 
to constant C 2 . 

Nevertheless, an informative assessment of constants Ci and C 2 
supplemented by an analysis of the the relationships 5 = tf(Re) re- 
sulted m another finding which simplified the whole task. By simple 
transformation, with a precision corresponding to the precision of 
the experiment, relation (15) can be transposed on Prandtl's known 
relation for a smooth tube. Such transformation corresponds to sub- 
stitution of the hydraulic diameter by a more suitable characteristic 
geometrical dimension. 

given pressure losses and flow-through weight of coolant ap- 
plies that 


5 = • rf/t 

Re = . dfo 
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If we designate a suitable characteristic dimension, which makes /27 
the mentioned transformation possible, as d 0 and the corresponsong va- 
lues of the coefficient of friction losses and the Reynolds criterion 
as £ 0 and Re 0 • then it applies that 


We denote 


€ o _ R&o _ do 

C Re. d\i 


(16) 


37 = H 

and call this value a geometrical factor. 
The 


(17) 


? o = (18) 

(19) 

and because the values C 0 and Re 0 apply to Prandtl's universal law of 
friction, i.e., 


I = 

/To" 

we can put it m the form of 


2 .log Re 0 


/Co - 0.0 


( 20 ) 


/H 

C 


2.log 


Re/C 


0.0 


(21) 


5.2 Heat Transfer 


Heat transfer into a channel with a variegated cross section can 
be described by means of a simple equation, such as, e.g., equation 
(3), only m certain cases and under certain conditions. Heat Trans- 
fer is generally described by equation 

q .1 • A f a U -01 dS (22) 


wherein /28 

<li = thermal output transferred from unit length 
of the heating element into the cooling gas 

•£ — local temperature of heat-exchanging surface 

0 = median temperature of gas 

■S = circumference of heat-exchanging surface 


We can simplify this equation on the basis of the analogy between the 
transfer of motion and heat. We assume that distribution t/T is, as 
well as distribution a/lx, merely a function of channel geometry. Ergo 


E = & (23) 

In regions with stabilized thermal conditions we can further assume 
that the prerequisite of automodel distribution of wall temperature 
has been met 


t - 0 , 

f 0-02 


LSI 


(24) 
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Equation (22) can then be written m the form 

q L = a. [t-e] (S) .,( 2 (S) .dS (25) 

The possibility for further simplification of equation (25) must 
always be considered in connection with the specific geometry of the 
channel and the specific design and material composition of the heated 
element. The justification for each simplification must be verified 
by preliminary calculations made under various assumptions and distri- 
butions of a along the heat-exchanging surface. Consideration in this 
case must be given to both internal as well as external factors (i.e., 
distribution of temperatures and thermal flows m the heated element) . 

The objective of these calculations is to find out in what way use of 
simplifications could distort energy calculations or the information 
which the experiment is to provide. 

Our subsequent comments involve only cases for which description 
of heat transfer could use an equation of the type 

q L = a(5>0j .S (26) /29 

where centering is done along the circumference of the heat-exchanging 
surface. 

Equation (3) used in assessment of thermokmetic experiments on 
geometry A and B represents application of equation (26) to our spec- 
ific geometrical conditions. The suitability for use of equation (3) 
as well as its limitations for channels with longitudinal fins was 
proven by a theoretical analysis and a number of calculations using 
the computer GIER. Details are contained in report [1]. 

Similarly as was the case with hydrodynamics, we first concentra- 
ted on the form of dependence of Nusselt's criterion on Reynold's cri- 
terion. We undertook a detailed analysis of our own as well as acces- 
sible experimental materiel of others and eliminated all cases in which 
we reached the conclusion that measurement was not performed at bal- 
anced hydrodynamic and thermal conditions. In a great majority of 
cases we then determined the dependence of 

Mu - contt. Re 0 * 6 * 8 (27) 

The constant is a function of channel geometry, thermophysical proper- 
ties of coolant and the thermal factor. If we further assume that the 
effects of thermophysical properties of the gasseous coolant can be 
expressed m the same way as for a tube, then relation (27) can be put 
in the form 

Mu = V . P/t° • 4 . Re 0 * 8 ..i|T n (28) 

where constant V is merely a function of geometry. 


6. Integral Geometrical Characteristics 


6.1 Pressure Losses 


The relation (21) generalizes our old experience that in certain 
cases involving noncircular channels the relation for a smooth tube /30 
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can be used for determination of £ provided that we use hydraulic dia- 
meter as a characteristic dimension of the channel, while m other in- 
stances this procedure could lead to misleading conclusions. 

Some authors tried to apply this experience m such a manner that 
m specific cases, inhere use of cfft did not facilitate unification of 
experimental results, they made use of another dimension of the given 
geometry as the characteristic dimension. As an example we can cite 
triangular lattices of rods with longitudinal bypass flow. If we use 
as a characteristic dimension the diameter of a circle plotted onto 
the lattice, the results show good coincidence with a smooth tube. 

Such an approach — if it produces results m the first place — is 
always dependent on the author's estimate and cannot be used without 
prior knowledge of experimental results. 

The kind of conclusions to which can lead an indiscriminate use 
of the hydraulic diameter and cntenal relation for a tube can be il- 
lustrated on three typical cases. 

Let us first consider a circular channel with very thin and deep 
internal radial slots. In such a case there is practically no differ- 
ence m the flow-through diameter in comparison with a smooth tube, 
while the wetted circumference can increase by a number of times. The 
hydraulic diameter and the friction loss coefficient will decrease m 
the same ratio. But it is clear that hydrodynamic conditions will not 
substantially differ from those in a smooth tube 1 . 

As a converse case we can take the eccentric annulus. According/31 
to the definition of dfr it should not be affected by accentricity and 
the latter should consequently also not affect pressure losses. How- 
ever, we know from the experiment that pressure loss decreases with 
increasing eccentricity and that the difference in pressure loss be- 
tween a concentric annulus and an annulus with eccentricity equal to 1 
amounts to 30 to 40%. However, hydraulic diameter in not an unequi- 
vocal characteristic value even for such simple geometry as is a con- 
centric annulus. Experiments have again proven that £ depends fairly 
heavily on curvature V\/Vz which is not reflected in dft. 

The relation (21) formally reduces this problem to the question 
of what determines the geometrical factor H and what does it depend 
on. That means that we must find a factor that would be merely a 
function of channel geometry and which would respect in a correct man- 
ner the indicated extreme cases. This factor must in some suitable 
way "weigh" the effects of individual parts of the wetted diameter on 
flow in the channel. If such a factor is to unequivocally define 
channels of the most varied geometrical configurations from the view- 
point of hydrodynamics, this factor must be of an integral character. 


1 For simplification's sake we have so far neglected general conditions 
that some authors (e.g., [13]) refer to as limiting the use of d( lt i.e., 

that in all parts of the flow-through diameter the width of the channel 
in the direction of the normal must be higher than the thickness of the 
"boundary" layer in which occurs maximum change in velocity during tur- 
bulent flow. 
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These general requirements could be used to define an infinite 
number of factors. We therefore embarked upon a search for a factor 
that would have the simplest and logical definition which would be, 
relatively speaking, the easiest for assessment and which could be 
applied at least in part to some descriptions of turbulent flux. 

After various considerations and carrying out a number of auxiliary 
calculations and unsuccessful experiments which we will not describe 
here, we decided on a factor that some authors use as a local measure 
of turbulence and which in the case of a circular channel is inter- 
preted essentially as the mixing length. 

The dissipative effect of an element of channel wall on random 
interference with flow in the vicinity of point M is quantitatively 
bound to reciprocal value of the distance of point M from the given 
wall element. For that reason the characteristic distance L of point 
M from the channel wall is given by the integral 

I'’! 2 { I d<f (29 > 

where 1 is the distance of point M from the channel wall in the direc-/32 
tion of cp. Relation (29) was used in his study [6] by Buleyev. Ac- 
cording this relation the factor L very closely resembles the measure 
of turbulence derived under assumption of resemblance of turbulent 
processes in various regions of liquid flux by A. M. Obukhov [40]. 

If we use relation (29) to assess the distribution of L in a circular 
channel we find out that 


J * 0.44 L 

The factor 1 in this case is the mixing path of Prandtl's semi- 
empincal theory of turbulent flow, arrived at by Nikuradze from ex- 
perimental values of velocity distribution in a circular channel with 
respecting the actual distribution of tangential stress in the tube 
cross section. 


The factor L, even though it respects the shape of the channel, 
is merely a local factor. However, the maximum value of L ( L max ) 
m the given cross section, or — what we deem more suitable — the 
median value of L in the given cross section could be taken as a 
value characterizing the channel. 

r « f p /L dV (30) 


This factor which has the dimension of width is not only a func- 
tion of channel configuration, but is at the same time a function of 
its absolute size. We therefore used as channel characteristic the 
factor 


which we shall refer to as 
acteristic. 



(31) 


the integral geometrical criterion or char- 


The factor 1 0 in relation (30) denotes the median characteristic 
distance T for an equivalent tube: 

T 0 = 0.08 87 . dfr (32) 
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6.2 Heat Transfer 


/33 


The flow-through cross section of noncircular channels can form 
an n-times coherent region bounded by n closed curves. Each of these 
curves corresponds to a part of the wetted circuit that can be indi- 
vidually heated. If m processing the values from investigation of 
pressure losses we could define channel geometry by a single geomet- 
rical characteristic, similar approach is not possible when dealing 
with heat transfer in a region with multiple connections. In such a 
case we must respect the distribution of a in the individual parts 
of the wetted circuit. (However, from the viewpoint of individual 
parts of the circuit, we again work with the median value of a m the 
given part — see equation (26)) . Thus, we need one more geometrical 
factor which would respect distribution of a. We undertook selection 
of this factor on the basis of the following simple contemplation: 


If we assume the validity of Reynold's analogy between transfer 
of motion and heat, then applies that 




where factor with the index "t" are related to the -c-th part of the 
wetted circuit. We further assume that no transfer in motion occurs 
accross the orthogonal trajectory to isotachs which intersects the 
relative maxima of velocity. This line is then plotted m approx- 
imation of the corresponding line in field L. 1 


If we now designate the area bounded by this trajectory and the /34 
heated circuit as Fj and the heated circuit Sj we can formulate 


a 



d 


. 4 Ft 

T ' ~sr 


We then used the factor dj/d k together with the factor L in process- 
ing experimental data regarding heat transfer m channels the flow- 
through cross section of which represents a doubly connected region. 


For a simply connected region it applies that 

d T =d h (34) 

£ 

Computation of both integral geometrical characteristics L and 
du/dj m the case of more complex geometries is practically impossible 
without a computer. Programs and calculations for our needs were pro- 
vided by Engr J. Schmid, CSc. Detailed methodology is m report [17]. 


1 We verified this assumption for the case of a smooth annulus, it being 
the only one for which experimental documentation is available. A part 
of the authors stipulates that the position of maximum velocity arrived 
at through assessment of the experiment conforms to the theoretical rel- 
ation for laminar flow also in the region of turbulent flow. Roberts in 
his theoretical study [18] uses for turbulent flow a relation derived 
from experiments, one that differs from laminar flow. The relative po- 
sition of maximum velocity m relation to the curvature of the annulus 
according to both relations is shown m fig. 26. The relative position 
of L ma>c lies between these two progressions. 
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7. Universal Crxtenal Relations 


I 

/35 j 

.1 

7.1 Universal Criterial Relation for C j 

1 

In assessment of the geometrical factor for the individual types 1 
of channels investigated we used the experimental relation ///£ = 

= *)(Re./tT) as a guide for determining whether the entire measurement or 
its part was affected by the transitory region £. Actual assessment of 
H was done graphically from experimental relations £ = ^ (Re) . With 
the use of a simple mechanical device (template) it was easy to per- 
form transformation to Prandtl's universally valid law of friction and 
directly read off the value of H. 

In our search for the dependence of H on channel geometry, the ex- 
perimental values of H were plotted consecutively in relation to L* , 
t-max/dj i and a combination of the two in the form ( L* nax /T . d^) . The re- 
sults are shown m fig. 21 through 24 in such a manner that. the indi- 
vidual investigated variants can be identified by markings of exper- 
imental points. (The used symbols are shown m tables for geometry 
A and B.) Figures 21 and 22 are practically identical, however the ■ 

latter figure shows only differentiation of geometry A and B and also 
includes results obtained by Rybacek [3]. j 

Smooth annulus is represented by a single point with indicated 
dispersion which includes the results obtained by the Institute of 
Nuclear Rese|rch, by Rybacek [3] and Ort [4], as shown in fig. 27. 

The factor L does not depend on it practically at all in the case 
of small annuli from a curvature of Vi/V 2 ~ 0.4 and up (see fig. 53) 
and has a value of 0.934. This value of L* represents the absolute 
minimum of L*. There is no geometry that could have a lower L * . 

< 

Comparison of fig. 21, 23 and 24 clearly shows that the relation 
W = f (L ) results in a substantially smaller dispersion of experi- 
mental points than is dependence of on two other factors. Thus, we j 

are of the opinion that the integral geometrical criterion L* defines 
channel geometry better than factors U max /d/ l ) and (Lf naK /I.df l ) . I 

The dependence of on channel geometry can be now expressed by 
the relation j 

H = 0.268 to 0.842 (L* ) “ 1 * 2 (35) 

Nevertheless, the exception in this relation is formed specif ical- /3( 
ly by a smooth tube, the law of turbulent friction of which we used as 
the basis for our processing. Relation (35) yield for a tube the value 
of H = 1.1 instead of the correct value H = 1. 

i 

We have failed so far in finding a satisfactory explanation for 
this contradiction. The simplest and most probable explanation might 
be that the used method, to include the integral geometrical character- 
istic, represents only a certain degree of approximation of actual con- 
ditions and, as such, cannot be in agreement with all empirical findings. 
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However, we must realize at the same time that a circular channel 
represents from the viewpoint of all other types of channels a certain 
extreme in two directions : 

1. It is the only channel in which the local value of tangential 
stress on the wall is identical with the median value of T. 

2. Of substantially more importance is the fact that for all 
practical purposes purely longitudinal flow exists solely in circular 
channels. In all other types of channels that were subjected to exper- 
imental investigation, there was evidence of a secondary transverse 
flow which strongly affects the distribution of t along the wetted cir- 
cuit and which necessarily must also become reflected in the extent of 
pressure losses. A possible exception from this allegation might be 
formed by annular channels. However, the results of measurements made 
so far still show enough divergence to preclude formation of an unequi- 
vocal conclusion. 

An answer to these problems could be provided only by a well-found- 
ed theory of turbulent flow m channels with a random cross section — 
if it can be formulated. 

Finally, plotted in fig. 25 are all our experimental results pro- 
cessed in the form of relation ( 21 ) with theoretical values of H ac- 
cording to relation (35) . The median error in comparison with rela- 
tion ( 20 ) is less than 10%, which in view of the overall accuracy of 
the experiments can be considered to be a very good result. 

In this manner the relation (21) in connection with relation (35) 
represents a criterial equation for the coefficient of friction losses 
the validity of which is not limited to merely a certain channel geo- 
metry. We will prove its potential application to all types of chan- 
nels sub para 8 . 


7.2 Universal Criterial Relation for Nusselt's Criterion /37 

All of our experiments were conducted in such a manner as to have 
the thermal factor in the region of balanced thermal conditions range 
between ^ = 1.2 to 1.3; median temperature ip = 1.27. As the character- 
istic thermokinetic factor for individual geometrical variants was se- 
lected and computed from experiments the value of nusselt's criterion 
for Re = 5.10 4 . We then gradually verified whether this value can be 
brought into an unequivocal dependence on some of the simple combina- 
tions of factors L* , L max and d( t /dj, which can be arrived on the basis 
of certain clear concepts of the transfer conditions m a channel. In 
this manner we have proven that this requirement can be met only by 
the geometrical criterion L* -A. The dependence of the values Mu .5 on 
the geometrical criterion L is shown in fig. 28. Marking of the 

experimental points again facilitates identification of the individual 
variants being investigated. 

Figure 29 shows the same dependence, but point markings serve only 
to differentiate geometry A and B. Experimental results arrived at by 
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K. Rybacek [3] are also included. His experiments were carried out 
at other values of ^ than were our experiments. As much as facili- 
tated by data m [3], _we computed the thermal conditions for the ex- 
periments done by Rybacek and converted his results to our conditions 
by means of Kutateladze ' s relation (13). In spite of the fact that 
the chosen approach can be considered to be only a rough approximation 
of actual conditions, the coincidence between those values and our own 
results is good. 

The interpolated relation in fig. 28 could be expressed by 

Mu = 97. - 91 (36) ' 

a T 

If we assume that even m the case of channels with a complex 
cmfiguration the effects of termophysical properties of the coolant 
can be expressed by means of the factor ?> i ° and the effect of the 
thermal factor by Kutateladze ' s relation, than from relation (36) we 
obtain a modified known relation for a circular channel: 

Mu = 0,011 . Pt 0 * 1 * * 4 * . Re 0 * 8 . ( L* 4)' 3 ’ 51 . tJ,- 9 * 5 5 (37) 

a T 

Thus, relation (36) precisely intersects the point which corres- /38 
ponds tothe value Mu 5 for an air-cooled circular channel. 

Let us return one more time to the diagrams in figures 28 and 29. 

In comparison with investigation of pressure losses, the results on- 
tained in experimental investigation of heat transfer with a view to 
the number and type of measured factors always burdened by larger er- 
rors. At the same time, investigation of heat transfer is more sub- 
jected to the distorting effects of the transitory region £ and of 
the entry sector. If we become aware of this fact, we can assess the 
dispersion of experimental points on these diagrams as relatively 
small. 


Figures 30 and 31 show plotting of all experimental relations 
Mu = £ (R&) adapted m the form 

Mu. (L*37)° 


. 9 1 


97 


= i (Re) 


(38) 


Points which lie beyond the region of ± 10% m comparison to the 
relation (37) belong for their major part to a few geometrical config- 
urations for which it was not possible to unequivocally determine 
whether the measurement was affected by one of the mantioned factors 
or not. 


Equation (37) can again be regarded as a critenal relation for 
Nusselt's number, which is not bound to a single specific geometry. 
Nevertheless, in comparison with the critenal relation for £, the use 
of this equation is limited by several conditions : 

1. The basic condition, which we pointed out already sub para 5, 

is suitability of equation (26) for computation of heat transfer in a 

specific case. Otherwise this condition may be defined so that heat 

transfer must be dependent on the actual distribution of the local co- 

efficient of heat transfer along the heat-exchanging surface. 


23 



2. Relation (37) was derived from experimental values measured 
in channels with a flow-through cross section represented by a singly 
of doubly connected region. No experimental documentation aws avail- 
able for multiple connected regions and suitability of relation (37) 
for such cases would first have to be established. 

3. Contrary to relation (21) , it was not possible to compare 
equation (37) with a larger number of other authors’ experimental 
findings regarding various channel geometries. 


7.3 Thermal Entry Sector / 

In an effort to make widest possible use of experimental materiel 
and to obtain as many specifically formulated findings, I attempted to 
derive from the experimental progressions of JL = ^ ( x ) the dependence 

of length Lj on channel geometry. We were interested in finding out 
whether integral geometrical characteristics could also be used in this 
case as an unequivocal geometrical criterion. We were aware at the 
same time of the fact that experimental documentation, in view of the 
overall concept of our experiments and accuracy of measurements, was 
not exactly the most suitable for this purpose. For that reason we 
regard our results to be merely the first bit of information offered 
about this problem and an incentive showing the approach to this par- 
ticular problem. Thus, we shall limit our comments to a mere listing 
of our conclusions: 


Distribution of the coefficient of heat transfer along the model 
can be described by equation 


- 2 - - / 


+ 0.1 


(-T i«i 


(39) 


This equation formally coincides with the equation for distribu- 
tion of a in a circular channel, which — with reference to the study 
of Boelter, Young and Iverson — is quoted m study [11]. 

The factor Lj ltl is an x-th section coordinate with = 1.1. 

This factor is a function of channel geometry and of Reynolds crite- 
rion, whereby dependency on Re is substantially weaker. In most vari- 
ants, the experimental results of which could be used for this purpose, 
*-T i.i increases with the Reynold's criterion up to values of Re = 2 to* 
3.10. For higher values of Re the value Lj 1 , remains constant. For 
these cases and for the region Re > 3.10 4 we devised the relations 

(40) 

* i (l *a£ <4i) 


The results are plotted in fig. 32 and 33. The geometrical char- 
acteristic L even m this case is not an unequivocal characteristic of 
the channel. As a result, relation (40) has one form for geometry A 
and another form for geometry B. 
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Relation (41) in comparison with the preceding relation is rela- /40 
tively unequivocal, but experimental points show considerable scatter- 
ln 9* From interpolation of experimental points we can express this 
relation by the equation 


j-T 1.1 
dh 



(42) 


For illustration we computed the experimental proqress of 
the form of relation 


a 

Otoo 


in 



• i ( 


Lt 


1.1 


The value Lt j j was for this purpose computed for individual var- 
iants from relation (42). The results are shown in fig. 34, which is 
being introduced without comment. 


7.4 Effects of Thermal Factor \p 

For some variants of geometry A we conducted experiments design- 
to P roduce results that would offer more accurate information about 
the effects of thermal factor <|> in channels with complex configuration. 

temperatures of the model were measured at a constant flow-through 
weight of coolant and at vaying values of thermal output and we deter- 
mined the coefficient of heat transfer in the region of developed 
thermal conditions. This approach was designed to suppress the effects 
of Reynold's criterion as much as possible. 

The results were processed in two ways: 

The first approach did not take into consideration any change 
m the Reynold's criterion in the investigated cross section as a re— 
suit of varying values of median gas temperature in individual regimes. 
Thus, effects m changes of Re on the Nusselt criterion were also dis- 
regarded. The results obtained through this approach are shown in fig. 
35 in the form of the relation 

Mu , , , . 

mnwr~n = 6 (rp} < 43 > 

To prevent an eventual distortion of the results as a consequence 
of an erroneous extrapolation of the relation M u = {(ip) into = i, W e 
related the relative value of Nusselt 's criterion to the value of Mu at 
y — 1.2. This value lies for all variants m the investigated range of 

r • 


.Et stands to reason that this approach is neither correct nor is /41 
it logical. Nevertheless, it is interesting that the experimental re- 
sults can be well described by Kutateladze' s relation (13). 

2. The second approach did take into consideration the effects 
of changes in gas temperature on the value of Re and the effects of 
the change on the value of W u. 
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The results are shown in fig. 36. Their analysis showed that 
for all the investigated variants the effects of the thermal factor 
can be expressed by the relation 


A lu 

AJu (4> = TJ 



(44) 


but that the exponent n is again a function of geometry. The depend- 
ence of n on channel geometry |.s expressed in two ways in fig. 37 and 
38.' One makes use of factor L , the other (L*^) as the geometrical 
criterion. 


In view of the conduct and precision of experiments and the used 
methodology for processing the results, point dispersion is not exces- 
sive in either case. However, there were entirely too few experiment- 
al points to warrant drawing of any conclusions from these results. 

As these experiments were conducted only on models of the geomet- 
ry Ai , it is not possible to judge which geometrical criterion would 
be more suitable for this purpose. However, it can be assumed that 
same as was the case with the thermal entry sector, the most unequi- 
vocal geometrical characteristic will turn out to be the criterion 

(*#■• 


8. Comparison of the Universal Criterial Relation for £ With Exper- /42 
imental Results Obtained by Other Authors 

To prove the universal applicability of the criterial relations 
(21) and (35) for the coefficient of friction losses m a wide range 
of channel geometries, we undertook a comparison of this relation with 
experimental values for all types of channels for which we could find 
credible documentation in technical literature. 


8.1 Triangular and Square Lattices of Parallel Rods 

A sheaf of parallel rods arranged into a triangular or square lat- 
tice may be regarded as one of the key geometries of channels from the 
viewpoint of nuclear power engineering. There is also a relatively 
large volume of experimental materiel available for this geometry. 

For these reasons, sheaves of rods with parallel by-pass flow received 
the greatest amount of our attention. 

Documentation for computation of H in both types pf an infinite 
lattice represents the dependence of the grometrical criterion L* or 
of the geometrical factor H on lattice density 6/d. The graphic form 
of these dependencies processed on computer GIER is shown in fig. 39. 
Comparison of experimental and computed values can be most vividly per- 
formed by comparing the relations 

S- - & (. 4 / d [ ( 45 ) 

wherein £o is the coefficient of friction losses in a circular channel 
at a given Re. The relation (45) can be substantially simplified for 
computation. The criterial relation for the coefficient of friction 
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losses (21) can be substituted in the region Re < 10 s by the relation 

| 0.376 (ijf)" 0 * 25 (46) 

After substitution from equation (46) into the relation (45) we obtain 

JL = H 1 25 (47) 

to 

wherein H = & (4/d). 

A comparison of this relation with experimental values for a tn- /43 
angular lattice is shown in fig. 41. It can be stated that there is a 
congruence of almost 100%. 

Investigations of velocity fields and of the distribution of tan- 
gential stress m an infinite triangulat lattice routinely make use of 
models with a configuration shown in fig. 42. However, due to their 
infinite dimensions these models are not suited for investigation of 
pressure losses, as has been experimentally shown a number of times 
already. Fig. 42 shows a comparison of experimental and computed 
values for this type of channel. For the sake of comparison, the same 
figure shows the theoretical dependence for an infinite lattice. The 
deviations of experimental values from theoretical dependence are 
smaller than is the dispersion of experimental points. Figures 43 and 
44 show a comparison of the experimental dependence of K = d ( R 2.) with 
theoretical dependence for various types of lattice density. Their 
coincidence is outstanding in all cases. Figure 44 also shows the 
potential extent of effects on overall losses that can be caused by 
inert draughts between sheaves of rods and the channel proper . 

The same conclusions, m regards to congruence between theoretic- 
al dependence and experimental results, can be drawn from diagrams 45, 

46 and 47, which show a comparison for the elements of a square lat- 
tice and an inifinite square lattice. 

Into this group of channels we could also include the channel in 
fig. 54. Three rods are placed in a triangular channel with rounded 
corners. Study [33] investigated the dependence of pressure losses 
in relation to the relative separation of rods. The forms of theoret- 
ical and experimental dependence are very similar in spite of the fact 
that dependence 5 * {\ (S/d) is not monotonous. However, with decreas- 
ing value of S/d the differences between the two progresses increase. 
Nevertheless, it is difficult to reach any conclusions from this com- 
parison. Experimental values were obtained for small values and a 
narrow range of Reynold's number. Figure 55 shows experimental re- 
sults recomputed into the relationship 7//? = & (Re/5). It is very 
difficult to determine from these progressions whether the measurement 
was affected by the transitory region of 5- However, that possibility 
should riot be excluded. 


8.2 Channels With a Triangular Cross Section 

Comparison of theoretical dependence with experimental values for 
channels with a cross section m the shape of an isosceles triangle 
turned out to go less in favor of theoretical dependence than was the /44 
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case in the preceding instance. Comparison of theory and experiment 
at varying values of the apex angle 3 of a channel is shown in fig. 

49 which has been essentially taken over from study [32]., In regard 
to form, theoretical and experimental dependence is mutually similar 
and at a value of 3 =60° they coincide also as to their absolute 
value. However, with decreasing 3 the difference between the two 
dependencies increases up to 12%. The progress of both dependences 
is interesting in the region 3 _i '0. We originally mistakenly assumed 
that this limiting case corresponds to parallel plates. Only compu- 
tation of the geometrical criterion L for a uniform trapezoid with 
varying values of the apex angle and varying ratio of height to base 
(fig. 52) showed how strongly changes the character of the dependence 
L = [t/Hj in the region 3*^0, t/k->0. 

The small divergence between theoretical and experimental pro- 
gress could be explained on the basis of the results obtained in the 
previously cited study [30]. The latter does not exclude the possi- 
bility that in the region of Re, in which the triangular channels 
were investigated, m a part of the channel could occur laminar flow. 
In that case the measurement results would be affected by the transi- 
tory region £. 

8.3 Sheaf of Rods With LongitudinalFins 

We undertook a comparison of the theoretical relation with ex- 
perimental results obtained by V. Krett [37] for a sheaf of rods, 
each with three longitudinal fins. The channel diagram and the mea- 
surement results are shown in fig. 56. Theoretical dependence is 
very well centered on experimental results obtained for isothermal 
and noni so thermal flow. 


8.4 Eccentric Annulus 

Eccentric annulus represents a typical example of a channel in 
which use of hydraulic diameter as a characteristic dimension fails. 

The more we were interested to find out how well would apply the cn- 
terial relation (21) with relation (35) . As we did not have at our 
disposal any studies providing experimental values, we used for com- 
parison the relation 

v-i— = 11 - e)° ’ 178 (48) 

^£ = 0 

(where is eccentricity) which was used in study [34] and was assess- /45 
ed on the basis of experimental findings in studies [35] and [36]. 

This relation does not take into consideration the experimentally 
detected weak effects of curvature. Its validity is also limited 
in the direction of an increase m e. 

In fig. 57 this relation is compared to two theoretical depend- 
ences for Vi /Vz - 0.5; 0.95. Congruence is again very good. 
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9. Conclusions 


/46 


1. An analysis of our and other authors' experimental findings 
in the measurement of friction and pressure losses in noncircular 
channels with hydraulically smooth walls showed that in all cases in 
which the results have not been affected by distorting factors, the 
dependence E, = £ (KeJ can be applied by simple transformation to the 
Prandtlian universal law of friction in a smooth tube. The transform- 
ation relation is hereby merely a function of channel geometry. 

2. On the basis of empirical findings we formulated general re- 
quirements for selection of geometrical characteristics factors which 
could play at a higher level a role similar to that played by hydraul- 
ic diameter at the outset of research into flux in noncircular chan- 
nels. On the basis of these requirements we defined integral geomet- 
rical characteristics (criteria) and it was empirically proven that 
these criteria relatively unequivocally define the geometry of non- 
circular channels from the viewpoint of pressure losses and transfer 
of heat. 


3. We used these geometrical criteria to process experimental 
materiel from research of pressure losses and of heat transfer in 
circular and annular channels with longitudinal fins. In this manner 
we obtained universal critenal relations for the coefficient of pres- 
sure losses and the Nusselt criterion. (Validity of the relation for 
W a is of course milited to gaseous media with Pa ~ 1.) 

4. Comparison of the criterial relation for E, with published ex- 
perimental values for noncircular channels showed applicability of 
this relation for channels with the most varied configuration with 
hydraulically smooth walls. 

5. Unification and generalization of experimental materiel which 
is facilitated by these criterial relations make it possible to expand 
our knowledge in regards to flux conditions in noncircular channels as 
regards integral characteristics and provide an improved basis for the 
study of turbulent flow theory than was offered previously by the wide- 
ly splintered experimental materiel. 
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Fig. 1. Geometrical conditions at entry into investigated channel - geometry A. 
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Fig. 2. 


Geometrical conditions at entry into investigated channel - geometry B 






Fig. 4. Dependence of the coefficient of friction losses in a smooth annulus on channel 
curvature according to various authors. (Taken over from study [4].) For key to this 
figure see table on next page. 
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Institute of Nuclear Research 
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Fig. 7. 

Dependence of 

C = i 

(Re) 

at varying conditions at entry 


into channels 

A2 and 

A3. 









Fig. 8. Dependence of £ = £ (Re) during varying conditions at entry 
into channels A10 and All. 
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Fig. 9. Dependence of f//f = i (Re/£) during varying 
conditions at entry into channels A2 and A3. 
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Fig. 10. Dependence of 7//£ = j$ (Re/£1 during varying 
conditions at entry into channels A10 and All. 







Fig. 13. Progress of a along the 
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Fig. 19. Several examples od experimental dependence of 
fju = i (Re) for geometry B. 
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Fig. 21. Assessment of the dependency of the geometrical 
factor on the geometrical criterion L from the experiments 
conducted by the Institute of Nuclear Research. 












Fig. 23. Dependence of geometrical factor on 
geometrical criterion (L max /d^). 
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Fig. 25. Processing of experimental results in the form of the relation 
£ _ / » Re. . with geometrical factor value acc. to equation (35) 
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comparison or maximum velocity in laminar 
and turbulent flow through smooth annuli 
with maximum value of L. 
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Values of the geometrical factor 
for varying curvatures of annular 
channels derived from experiments, 
and progress of H = &[V 1 /t> 2 ) ac- 
cording to relation (35). 
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Fig. 29. Experimental values of Nusselt’s criterion for Re = 5.10f* and if = 1.27 
and their dependence on the geometrical criterion (L } differentiation between 
geometry A and B, including results obtained by K. Ryfiadek [3]. 
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Fig. 30. Processing of the results of thermokinetic experiments on geometry A into the 
dependence form _!_{ fjk) °* 91 . Uu = {, (Re)- 
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Fig. 31 Processing the results of thermokinetic experiments on geometry B into the 
dependence form _!_( L * # Wa = ^ ( 4e ) . 
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Fia 34- Distribution of a along the length of the model m the form of 
dependence^- - 1 = £ (-q^-) . Values of L Tm for individual 

variants were computed ^rom equation (42). 
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fig. 38. Experimental values of exponent n as a function of criterion L. 
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Fig. 42. Comparison of computed and experimental 
values of 5/5o for the used cell models of an infi- 
nite triangular lattice. 
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Fig. 43. Comparison of experiments and computations for a triangular 

when it can be considered to be infinite) . 
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Fig. 44. Continuation of Fig. 43. 









































Fig. 48. Dependence of L and H on the apex angle of an isosceles triangle 



Fig. 49. Comparison of computed and experimental values [32] for a 
channel with a cross section in the shape of an isoscles 
triangle. 
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Fig. 50. Coefficient of friction losses 
m an isosceles triangle [29]. 




Fig. 51. 


Extent of the laminar region m a channel with an 
triangle cross section with an apex angle of 11.5° 


isosceles 

[30]. 
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Fig. 54. Comparison of computed and experimental data of study [33], 
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Fig. 55. Results of study [33] recomputed into form ///£ 
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GEOMETRICAL CHARACTERISTICS Of MODELS 


Variant // No .. f 
& symbol ot 

fins 

ttnocent. n 


TABLE I. 


— 


Main . Outer dia of inner dia of Height 


channel modol model tube 


of fin 

h 


3 

H 

4 

+~ 

5 

X 

6 O 

7 

□ 

8 

O 

9 

A 

10 , 

A 

n . 

A 

12 

A 

13 

Y 

14 

A 


81, 0 + 0,1 
M 


70,3 + 0,04 


8 qq + 0 * ^ ^ 

3 41 * - 0 ,?5 




>r 

19 

41 

41 ’ -0,03 


0,3 

0,25 

o,T 

0,02 

b,on 

0, '2 

0,16 

0,14 

0,25 

0,-°5 

0,03 

0,07 

0,2 

0,2 

o'.TJ 


0 41 
3 41 


I 29 + ^ 

II -0,09 



































































































GEOMETRICAL CHARACTERISTICS OF MODELS B 



4 


' m. -• 

! • +* • 

m 


3 3 3 fc £ IT 
* ~ □ * 0 




B • ID 

i 


a >.> © 

o 




© fc- B ® > © 

R £ 3 X 3? F 




Key: tl) Variant #, & symbol; (2) Number of fins; (3) Inner diameter; 

(4) Outer diameter; (5) Height of fins; (6) Tickness of fins; (7) Flow- 
through area; (8) Wetted circuit; (9) Remark; (10) Not investigated; 

(11) Shape of fins considerably irregular. 



Table of Annular Channels 
Measured at the Institute of Nuclear Research 


Designation ^2 


Ml 


00,98 


80,98 


70,? 


70,3 


63,99 


63,99 




50,65 


60,55 


50,65 


60,55 


50,65 


60,55 


V 10 W 


3 ,o33 


?,043 


1 ,965 


0,975 


1,334 


0,344 


0,625 


0,740 


0,7? 


0,861 


0,792 


0,946 


Geometrical Characteristics of Models from study [3] 
SMOOTH ANNULAR CHANNELS 


Type of cooling V d T 


I>„ 


am mm 


of model 


2,376 

1,52 


38 32 1,183 2,0 


6 2,176 2,0 


?? 0,93? cooled core 

13 0,932 

6 0,935 


22 0,932 bilateral cooling- 


38 16 2,376 2,0 


22 0,932 cooled sleeve 1,06 


ANNULAR CHANNELS WITH A LONGITUDINALLY FINNED CORE (GEOMETRY A) 


Length 


mm mm 


45,3 


12 45,8 


of model 


ra 


1,996 


6,22 

1,26 

5,73 

1,05 

6,22 

1,?6 


15,75 


12,05 


10,52 

6,81 


6,79 


1,108 


1,000 


1,036 


1,019 


1,106 


1 

1 

0 

0 
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